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Abstract 

■ In this paper we explicitly work out the secular perturbations in- 
duced on all the Keplerian orbital elements of a test body to order 
C(e 2 ) in the eccentricity e by the weak-field long-range modifications of 
the usual Newton-Einstein gravity due to the Dvali-Gabadadze-Porrati 
(DGP) braneworld model. Both the Gauss and the Lagrange pertur- 

| bative schemes are used. It turns out that the argument of pericentre 

t^J- ■ u> and the mean anomaly M. are affected by secular rates which depend 

' on the orbital eccentricity via 0(e 2 ) terms, but are independent of the 

in , 

semimajor axis of the orbit of the test particle. For circular orbits the 
Lue-Starkman (LS) effect on the pericentre is obtained. Some observa- 
tional consequences are discussed for the Solar System planetary mean 
longitudes A which would undergo a 1.2 • 10~ 3 arcseconds per century 
5h i braneworld secular precession. According to recent data analysis over 

J^P| 92 years for the EPM2004 ephemerides, the 1-sigma formal accuracy in 

^ i determining the Martian mean longitude amounts to 3 • 10 -3 milliarc- 

seconds, while the braneworld effect over the same time span would 
be 1.159 milliarcseconds. The major limiting factor is the 2.6 • 10 -3 
C$ ' arcseconds per century systematic error due to the mismodelling in 

the Keplerian mean motion of Mars. A suitable linear combination of 
the mean longitudes of Mars and Venus may overcome this problem. 
The formal, 1-sigma obtainable observational accuracy would be ~ 7%. 
The systematic error due to the present-day uncertainties in the solar 
quadrupole mass moment J2, the Keplerian mean motions, the general 
relativistic Schwarzschild field and the asteroid ring would amount to 
some tens of percent. 
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1 Introduction 

Recently, a braneworld scenario which yields, among other things, long- 
range modifications of the Newton-Einstein gravity has been put forth by 
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Dvali, Gabadadze and Porrati (DGP) ^H]- In such model, which encom- 
passes an extra flat spatial dimension, there is a free crossover parameter 
ro, fixed by observations to a value ~ 5 Gpc, beyond which traditional grav- 
ity suffers strong modifications yielding to cosmological consequences which 
would yield an alternative to the dark energy in order to explain the ob- 
served acceleration of the Universe. The consistency and stability of the 
DGP model have recently been discussed in 

Interestingly, for R g « r « ro, where R g = 2GM/c 2 is the usual 
Schwarzschild radius for a gravitating body of mass M, there are also 
small modifications to the usual Newton-Einstein gravity which could be 
detectable in a near future. Indeed, Lue and Starkman (LS) derived in 
[I] an extra-pericentre advance for the orbital motion of a test particle as- 
sumed to be nearly circular. Its magnitude is ~ 4 • 10 -4 arcseconds per 
century (" cy _1 in the following). In |3J it has been shown that the sign 
of such an effect is related to the cosmological expansion phases allowed 
in this model: the Friedmann-Lemaitre-Robertson- Walker (FLRW) phase 
and the self-accelerating phase. The LS precession is an universal feature of 
the orbital dynamics of a test particle because, in this approximation, it is 
independent of its orbital parameters. 

Since the present-day accuracy in measuring the non-Newtonian perihe- 
lion rate of Mars amounts to ~ 10~ 4 " cy _1 (E.V. Pitjeva, private communi- 
cation 2004), the possibility of measuring such an effect in the Solar System 
scenario from the planetary motion data analysis seems to be very appeal- 
ing. It has been investigated with some details in iij] where a suitable linear 
combination of the perihelia of some inner planets has been considered. 

In this paper we work out the secular effects of the DGP model on all the 
Keplerian orbital elements of a test body to order 0(e 2 ) in the eccentricity e 
with the Gauss and Lagrange perturbative schemes. Possible observational 
implications are worked out. 

2 The orbital effects 

In this Section we will work out the secular effects of the DGP gravity on 
the Keplerian orbital elements of the orbit of a test body which is depicted 
in Figure ^ 

From the metric for a static, spherical source in a cosmological de Sitter 
background @] 

(ds) 2 = A^ 2 (r,u;)(c(it) 2 -^ 2 (r,u>)(dr) 2 -5 2 (r,w)[(^) 2 +sin 2 0(#) 2 ]-(du;) 2 , 

(1) 
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Figure 1: Orbital geometry for a motion around a central mass. Here L 
denotes the orbital angular momentum of the particle of mass m, J is the 
proper angular momentum of the central mass M, II denotes the pericentre 
position, / is the true anomaly of m, which is counted from II, f2, to and 
i are the longitude of the ascending node, the argument of pericentre and 
the inclination of the orbit with respect to the inertial frame {x, y, z} and 
the azimuthal angle <f> is the right ascension counted from the x axis. When 
orbits with small inclinations are considered, the longitude of pericentre 
w = ft + ijj is used. 
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where w is the fourth spatial coordinate, and 



N ~ 1 



2r 



the following Lagrangian can be obtained for a spherically symmetric matter 
source located on the brane around the origin 1 (r = w = 0) 



£dgp 



m 
~2 




(*) 2 - (y) 2 - (i) 2 



(3) 



2.1 The Gauss perturbative scheme 

From 

d_ (d£\ _ dC 
dt \dv J dr 

the resulting braneworld acceleration is 



iDGP 



GM 



(4) 



(5) 



where r is the unit vector in the radial direction. The term © can be 
regarded as a small perturbation whose effects on the Keplerian orbital 
elements of a test particle can be straightforwardly worked out, e.g., in 
the Gauss perturbative scheme (e.g. |Hj). Note that © is purely radial. 

The Gauss rate equations for the semimajor axis a, the eccentricity e, 
the inclination i, the longitude of the ascending node O,, the argument of 

1 Since we are in the weak-field approximation, it is assumed N(r,w) = 1 + 
h(r,w), A(r,w) = l + a(r,w), B(r, w) = r[l + b(r, w)]. The gauge b(t, r)\ w= o = has been 
adopted. The correction only to the Newtonian potential has been considered, i.e. the 
spatial part of the metric has been considered as Euclidean. Note that the modifications 
in the cosmological background introduced in do not alter these results. 
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pericentre u> and the mean anomaly M of a test particle are 
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in which n = 2tt/P = sj GM/a 3 is the Keplerian mean motion, P is the test 
particle's orbital period, / is the true anomaly counted from the pericentre, 
p = a(l — e 2 ) is the semilactus rectum of the Keplerian ellipse, Ar, At, An 
are the in-plane radial, transverse and the out-of-plane components of the 
perturbing acceleration, respectively, which have to be evaluated on the 
unperturbed Keplerian ellipse 



a(l 



1 + e cos / 



(12) 



The secular effects can be obtained by averaging over one orbital period the 
right-hand-sides of ((S|l-(|11|) evaluated on (fT2*)) . The average is performed by 
using 



dt 
P 



1 



l f/ 2 df 



2vr(l + ecos/) 2 (13) 

and integrating from to 2ir. Since the perturbing acceleration is entirely 
radial, from ©-© it is straightforward to note that there are no perturba- 
tions in the inclination and the node. In regard to the other elements, the 
following expansions have been performed 

(1 + ecos/r 3 / 2 ~ 1- 3 ecos/ + -e 2 cos 2 /--e 3 cos 3 /, (14) 

2 8 16 



(1 + ecos /) 



-5/2 



i " t i ^ 2 2 t 

1 ecos t H e cos t. 

2 J 8 J 



(15) 
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The expansion (|14[) has been used for the semimajor axis, the eccentricity 
and the pericentre, while the expansion (|15|) has been used for the mean 
anomaly. 

It turns out that only the pericentre and the mean anomaly undergo 
secular effects which are, to order 0(e 2 ) 

du 3c ( _ 13 2 



dM , lie A 39 o 



±— 1 - — e l . (17) 
dt 8r V 352 y v 7 

For e — > lfTH|l reduces to the result by Lue and Starkman for circular orbits 

EB- 

2.2 The Lagrangian perturbative scheme 

The Lagrange planetary equations for the rates of change of the Keplerian 
orbital elements are (e.g. jUJ) 



da 2 dTZ 



dt na dM ' 

de (1 - e 2 ) dTZ VT^e 2 dTZ 



(18) 
(19) 



dt na 2 e dM na 2 e duo ' 

— cosi 1 — 1 dTZ ^ 

dt na 2 Vl — e 2 sini <9w na 2 Vl — e 2 sini 917 

do i m (21) 



rea 2 Vl — e 2 sini <9z 

i an VT^e 2 an . . 

— = -COSI ; + - o ^> 22 

dt na 2 Vl-e 2 sinz di na z e de 

— (i-e 2 ) 971 2 <^ 

na 2 e de na da ' 

where TZ is the perturbing function which accounts for all the departures 
of the gravitational potential from the Newtonian monopole term. Such a 
perturbing scheme was applied for the first time to post-Newtonian motions 
by Rubincam ^U]. In the case of the DGP braneworld theory from Ti = 
p • v — C, p = dC/dv, and © it can be obtained 

K DGP = =F— V GMr. (24) 
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The DGP perturbing function (|24(1 has to be evaluated on the unperturbed 
Keplerian ellipse (|12j) and averaged over one orbital revolution by means of 
(|13f) in order to obtain the secular effects. From 

(1 + e cos f)~ 5/2 ~ 1 - -e cos / + — e 2 cos 2 / - — e 3 cos 3 / + H^ e 4 cos 4 f 
v J ' 2 J 8 16 128 

(25) 

one obtains 

cVGMa A , 3 2 9 4 



w ~^^l 1 + n e + i^ e J- <26) 

Since 

is it straightforward to obtain 

da de di dQ , . 

— = — = — = — = 0. 28 
dt dt dt dt v ' 

The situation is different for the pericentre and the mean anomaly. Indeed, 



d{K) 3ecVGMa~ ( 3 2 



d(K) c GM ( 3 2 9 , 



da 2r V a \ 16 1024 

As a consequence, 

cftj 3c A 13 2 



1 + TS e + T^7 e ■ (30) 



dM , 11c A 39 2 



±_ i _ . (32) 

8r V 352 y v ; 

Note that (f5T|)-(|52|) are identical to (fTB |) -lfT7 |l obtained in Section |2~T1 



3 The possible use of the planetary mean longi- 
tudes 

In [H] only the use of the planetary perihelia has been examined. First, 
preliminary observational tests are discussed in [7]. 

The fact that also the mean anomaly is affected by the DGP gravity 
with a relatively large effect suggests to examine the possibility of using the 
planetary mean longitudes A = to + $7 + M. . 
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3.1 The observational sensitivity 

The sizes of the DGP secular precessions and of many Newtonian and Ein- 
steinian competing secular precessions of the mean longitudes of the Solar 
System planets are reported in Tabled It can be noted that the DGP secu- 
lar precessions amount to 1.2- 10~ 3 " cy _1 for all the planets: the corrections 
due to the eccentricities are very small amounting to almost 10~ 5 " cy _1 . 
In regard to the observational sensitivity, in Table El retrieved from Table 4 
of the most recent results for the EPM2004 ephemerides are presented. 
They are based on the processing of a vast amount of data of different kinds 
(optical, radar, spacecraft, LLR) ranging from 1911 to 2003. It can be noted 
that Mars is the best candidate for extracting the DGP effect because the 
formal standard deviation in AMars amounts to 3- 10~ 3 milliarcseconds (mas) 
only while the braneworld shift for the same time span is 1.159 mas. 

4 Some systematic errors 

In this Section we will examine various competing classical and general rel- 
ativistic effects which would act as sources of systematic errors in order to 
see if it is possible to use only the Martian mean longitude for the proposed 
test. 

4.1 The impact of the solar quadrupole mass moment and 
of the Einstein gravitoelectric force 

In regard to the systematic errors which would induced by the other Newto- 
nian and Einsteinian competing effects, the solar quadrupole mass moment 
J2 is presently known with a 15% accuracy (Table 7 of 11 ), so that the 
mismodelled part of its secular precession would be two orders of magnitude 
smaller than the effect of interest. 

More important would be the impact of the Einsteinian gravitoelectric 
precession. Indeed, the general relativistic perihelion precession has been 
measured to a 10 -4 relative accuracy via the PPN parameters (3 and 7 
(Table 8 of ^J): assuming that it would also hold for A, this would yield a 
2 • 10~ 4 " cy _1 mismodelled effect. However, it should be noted that these 
are merely the formal 1-sigma errors: realistic bounds might be one order 
of magnitude larger. 
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Table 1: Nominal values, in " cy~ , of the secular precessions in- 
duced on the planetary mean longitudes A by the DGP gravity and by 
some of the competing Newtonian and Einsteinian gravitational pertur- 
bations. For a given planet, the precession labelled with Numerical in- 
cludes all the numerically integrated perturbing effects of the dynam- 
ical force models used at JPL for the DE200 ephemerides. E.g., it 
also comprises the classical N- body interactions, including the Keple- 
rian mean motion n. For the numerically integrated planetary precessions 
see http://ssd.jpl.nasa.gOv/elem_planets.html#rates. The effect labelled 
with GE is due to the post-Newtonian general relativistic gravitoelectric 
Schwarzschild component of the solar gravitational field, that labelled with 
J2 is due to the classical effect of the Sun's quadrupole mass moment J2 
and that labelled with LT is due to the post-Newtonian general relativistic 
gravitomagnetic Lense-Thirring ^21 component of the solar gravitational 
field (not included in the force models adopted by JPL). For J2 the value 
1.9 • 10" 7 has been adopted For the Sun's proper angular momentum 
J, which is the source of the gravitomagnetic field, the value 1.9 • 10 kg 



m 2 s 1 __3I has been adopted. 


Planet 


DGP 


Numerical 


GE 




J 2 


LT 






Mercury 


1.2 • 10" 3 


5.381016282 • 10 8 


-8.48 • 


10 1 


4.7- 10" 2 


-2 


10' 


-3 


Venus 


1.2 • 10~ 3 


2.106641360 • 10 8 


-1.72 • 


10 1 


5 • 10~ 3 


-3 


10' 


-4 


Earth 


1.2 • 10~ 3 


1.295977406 • 10 8 


-7.6 




1.6- 10~ 3 


-1 


10" 


-4 


Mars 


1.2 • 10~ 3 


6.89051037 • 10 7 


-2.6 




3 • 10~ 4 


-3 


10- 


-5 


Jupiter 


1.2 • 10~ 3 


1.09250783 • 10 7 


-1 • 10 


-1 


5 • 10~ 6 


-7 


10- 


-7 


Saturn 


1.2 • 10~ 3 


4.4010529 • 10 6 


-2- 10 


-2 


6 • 10~ 7 


-1 


10" 


-7 


Uranus 


1.2 • 10~ 3 


1.5425477 • 10 6 


-4- 10 


-3 


5 • 10~ 8 


-1 


10- 


-8 


Neptune 


1.2 • 10~ 3 


7.864492 • 10 5 


-1 • 10 


-3 


1 • 10- 8 


-5 


10- 


-9 



Table 2: Formal standard deviations, in milliarcseconds (mas), of the plane- 
tary mean longitudes as from Table 4 of for the EPM2004 ephemerides. 
Note that realistic errors may be an order of magnitude larger. About 
300000 position observations (1911-2003) of different types (optical, radar, 
spacecraft, LLR) have been used. The braneworld shift for A amounts to 
1.159 mas over a 92-years time span. 



Mercury 


Venus 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


3.75 • 1Q- 1 


1.87- 1Q- 1 


3 • 10~ 3 


1.109 


3.474 


8.818 


3.5163 • 10 1 
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4.2 The N-body perturbations 

The mean longitude can be written as A = nt + e, so that ^5] 

dX ^ de* _ ^ dn ^ ^ de 
dt dt dt dt ' 

where 



de _ 2 yT^g(l - yT^g) dK | tan(i/2) d7£ 

na da na 2 e de na 2 ^/l — e 2 <9i 

4.2.1 The direct and indirect effects on the Keplerian mean mo- 
tion due to the other planets and the asteroids 

In regard to the Keplerian mean motion, it turns out to be the major limiting 
factor. Indeed, the uncertainty in the solar GM ugMq = 8 • 10 9 m 3 s -2 
(Table 1 of Q3|) yields 

""Mars — 2-6 -10 3 " cy . However, if we adopt 
the position of keeping fixed the solar GM, the formal 1-sigma error due to 
the semimajor axis cr aMara = 6.57 • lO^ 1 m (Table 4 of only reduces to 
^ Mafs =3-10- 4 "cy- 1 ! rS 

Another source of potential bias is represented by the indirect pertur- 
bations on the Keplerian mean motion induced by the variations in the 
semimajor axis 



dn 3nda 3 d(TZ) 

An = —t = -t= — o-^t. 35 

dt 2adt a 2 d\ 



According to 15 , there are no secular effects on the semimajor axis: instead, 
the so called resonant perturbations, affect this Keplerian orbital element. 
They are induced by those terms in the expansion of the N-body disturbing 
function which retain the mean longitudes of the perturbed and the perturb- 
ing bodies. Such kind of harmonic perturbations, due to the asteroids for 
Mars, are potentially very insidious because they may have large amplitudes 
and extremely long periods. This topic has been treated in ^Hj. Table IV 
of j!6j lists the most important of such perturbations. The nominal ampli- 
tudes of the perturbations induced, e.g., by (1) Ceres, (2) Pallas, (4) Vesta 
and (7) Iris are 4.7 • 1(T 2 " cy -1 , 1.2 • 10~ 2 " cy" 1 , 5.7 • 10~ 2 " cy" 1 and 
5 • 1CP 3 " cy -1 , respectively. According to the results of Table 6 of m, their 
mismodelled parts amount to, 7- 10 -5 " cy -1 , 3- 10 -5 " cy -1 , 4- 10 -5 " cy -1 
and 8 • 10~ 5 " cy -1 , respectively. It must also be noted that the integrated 
shift of An grows quadratically in time. 



10 



4.2.2 The perturbations on the mean longitude due to the other 
planets and the asteroid ring 

Here we will deal with de/dt which is responsible for the secular perturba- 
tions on A. 

The Newtonian secular perturbations induced on the Mars mean lon- 
gitude by the other planets of the Solar System are of the order of 10 3 " 
cy -1 . The major source of uncertainty is represented by the GM of the 
perturbing bodies among which Jupiter plays a dominant role, especially 
for Mars. According to |17) . the Jovian GM is known with a relative ac- 
curacy of 10 -8 ; this would imply for the red planet a mismodelled pre- 
cession induced by Jupiter of the order of 10~ 5 " cy^ 1 . The GM of Sat- 
urn is known with a 1 • 10 -6 relative accuracy ^E]- However, the ratio of 
the secular precession induced on the Martian Keplerian elements by Sat- 
urn to that induced on the Mars perihelion by Jupiter is proportional to 
(^Sat/Mjup)(ajup/asat) 3 ~ 5 • 10~ 2 . This would assure that also the effect 
of Saturn is of the order of 10 -5 " cy . The situation with the precessions 
induced by Uranus and Neptune is even more favorable. Indeed, for Uranus 
(Mura/Mjup)(aj U p/aura) 3 ~ 1 • 10~ 3 and the relative uncertainty in the ura- 
nian GM is 2 • 10" 6 |Ej|. For Neptune (M Ncp /Mj up )(aj up /a Nep )~3 • 10 -4 
and ctgm/GM = 2 • 10" 6 [10]. 

A source on potentially non-negligible perturbations on the Martian 
mean longitude is the asteroid ring, i.e. the ensemble of the minor as- 
teroids whose impact can be modelled as due to a solid ring in the ecliptic 
plane The perturbations due to it can be worked out, e.g., with the 

Lagrangian approach and the disturbing function of the Appendix of |21| . 
By using the values of for the ring's radius and mass it turns out that 
the secular perturbation on AMars amounts to —3.4 • 10 -3 " cy -1 , with an 
uncertainty of 3 • 10 " cy -1 . 

4.3 The total systematic error on the mean longitude of 
Mars 

In Table El we summarize the various systematic errors affecting the mean 
longitude of Mars. It turns out that it is impossible to only analyze the 
mean longitude of Mars: indeed, the 1-sigma total error would amounts to 
~ 1.1 • 10~ 3 " cy j i.e. more than 90% of the braneworld effect. 

In the next Section we will outline a possible strategy to suitably combine 
the data of Mars with those of other inner planets in order to cancel out, 
by construction, many systematical error. A numerical example is explicitly 
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Table 3: Sources of systematic errors, in " cy -1 , affecting the mean longitude 
of Mars. The reported figures are at 1-sigma level. The total effect, obtained 
by summing up the various errors, is more than 90% of the braneworld 
signature. 

Source of systematic error Mismodelled amplitude (" cy -1 ) 

Keplerian mean motion 3 • 1CP 4 

Asteroid ring 3 • 1CP 4 

Schwarzschild GE field 2 • 1(T 4 

(7) Iris 8 • 1(T 5 

(1) Ceres 7 • 1(T 5 
Solar J 2 4.5 • 1(T 5 
(4) Vesta 4 • 1(T 5 

(2) Pallas 3 • 1(T 5 
Lense-Thirring GM field (assumed unmodelled) 3 • 10~ 5 

Total 1.1 • 1(T 3 



worked out. 



5 The linear combination approach 

In order to cancel out the impact of the various sources of systematic errors 
it is possible to suitably linearly combine the mean longitudes of Mars and 
Venus following an approach adopted in, e.g., |22l|2Sj. Let us assume we 
have at our disposal N Keplerian orbital elements 2 K, whose time evolution 
is supposed to be affected by a certain number of Newtonian and post- 
Newtonian effects, say 

K- = ^DGP + K-,h + ^GE + ^N-body + K-LT + ■■■■ (36) 

If we are interested in isolating one particular feature, say Adgp> and we 
know it is smaller than other larger effects which affect the same Keplerian 
element we can explicitly write down the expressions of the observational 
residuals 3 5IC Q ^ S in term of the feature of interest- which will be assumed 
to be entirely (or partly) present in the residuals-and of the main larger 

2 They could all belong to the same planet or, alternatively, they could all be the same 
element, say, the mean longitude, of N planets or a mix of these possibilities. 

3 Here we speak about residuals of Keplerian orbital elements in a, strictly speaking, 
improper sense. The Keplerian orbital elements are not directly observable: they can only 
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aliasing effects-which will affect the residuals with their mismodelled part 
only-so that the number of terms in the sum in the right-hand-side of (|36jl 
which represent the effect of interest and the other most relevant larger 
bias is equal to the number N of Keplerian orbital elements we have at our 
disposal. Now we have a system of N equations in N unknowns which we 
can solve for the effect we are interested in. The resulting expression will 
be, by construction, independent of the other larger aliasing effects. 

By using, e.g., the figures of Table [Qhic the numerically integrated total 
precessions, which encompasses all the Newtonian and general relativistic 
effects, it is possible to obtain 

<5A Mar s + ci^Avenus = ^8 ■ 10~ 4 "cy~\ (37) 

where 

; (Num) 

cx = Jjs™ = _3.270 • 1CT 1 

; (Num) 
A Mars 

and 5X are the time series residuals of the mean longitude built up in order 
to entirely absorb all the non-Newtonian and non-Einsteinian gravity. 

The combination (|37|) is not affected, by construction, by all the com- 
peting Newtonian and general relativistic perturbations acting upon A, at 
least to the level of accuracy of the dynamical force models used in calculat- 
ing A( Num ) . Of course, the coefficient c\ can, in principle, be more precisely 
recalculated by using future, more accurate ephemerides. Indeed, from the 
previous discussion should be clear the importance of also including, at the 
best of our knowledge, the asteroids. While in the mathematical model of 
DE200 the perturbations of only three of the major asteroids are present, in 
the more advanced DE410 and EPM2004 ephemerides the perturbations of 
300 asteroids and also of the asteroid ring have been included. 

From Table I2 it is possible to obtain for the formal 1-sigma observational 
accuracy in (|37|) a 7% value over 92 years. The systematic error due to the 
solar quadrupole mass moment, the Keplerian mean motions, the general 
relativistic Schwarzschild field and the asteroid ring amounts to some tens 
percent by assuming the formal, 1-sigma level of uncertainty of for J2, 

AMars, O Venus, P, 7 an d Af ring . 

be computed. The basic observable quantities are ranges, range-rates and angles. Here we 
mean the differences between the time series of K. obtained from a given observed orbital 
arc and the time series of K, obtained from a propagated orbital arc with the force we 
want to cancel in the force models. The two time series share the same (measured) initial 
conditions. 
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(38) 



Table 4: Orbital parameters of LAGEOS, LAGEOS II and LARES. 



Orbital element 


LAGEOS 


LAGEOS II 


LARES 


semimajor axis a (km) 


1.2270 • 10 4 


1.2163 • 10 4 


1.2270 • 10 4 


eccentricity e 


4.5 • 10~ 3 


1.4 • 10~ 2 


4.0 • 10" 2 


inclination i (deg) 


1.10 • 10 2 


5.265 • 10 1 


7.0 • 10 1 



6 The possibility of using the proposed LARES /WEBER- 
SAT satellite 

The Earth artificial satellite LARES/WEBER-SAT PHH] was proposed in 
order to measure the Lense-Thirring effect on the orbit of a test particle to a 
high accuracy level (~ 1%) in the gravitational field of the Earth by suitably 
combining its data with those of the existing geodetic SLR (Satellite Laser 
Ranging) satellites LAGEOS and LAGEOS II. Their orbital parameters are 
in Table 0J In [2f)j Ciufolini has proposed to measure the braneworld effect 
on the pericentre, which amounts to ~ 4- 10 -3 milliarcseconds per year (mas 
yr~ 4 ), with the perigee of LARES. Unfortunately, it would be impossible, 
as we will show in the following. 

6.1 The observational sensitivity 

From Ar ~ eaAu; [57j, it can be obtained that the accuracy in measuring 
the perigee precession over a given observational time span can be expressed 
as Su ~ 5r/ea. By assuming a root-mean-square (rms) error of 1 mm 
in reconstructing the LARES orbit over, say, one year for a given set of 
dynamical force models one gets 5uj = 4 • 10 _1 mas. 

With 4 a = 3.6 • 10 4 km and e = 2.8 • 10 _1 , as also suggested in [2Hj, and 
by assuming a rms error of 1 cm, the accuracy in the perigee would amount 
to 2 • 10 _1 mas. Note that the adopted values for the obtainable accuracies 
in r are optimistic; for example, the mm accuracy has not yet been fully 
achieved for the existing LAGEOS satellites. 

6.2 The systematic errors of gravitational origin 

The perigee of an Earth artificial satellite is affected by various kinds of 
long-period (i.e. averaged over one orbital revolution) orbital perturbations 

4 It may be interesting to note that, with such an orbital configuration and i = 6.34- 10 1 
deg the gravitomagnetic Lense-Thirring precession would amount to -1.8 mas yr only. 



14 



induced by the multipolar expansion of the Earth's gravitational potential 
[S]. The most insidious ones are the secular rates induced by the even 
(£ = 2,4,6...) zonal (m = 0) harmonic coefficients Jg of the geopotential 
which account for the departure of the Earth from an exact spherical shape. 
Their explicit expressions up to degree I = 20 can be found, e.g., in 28 . 
The largest precession is induced by the Earth's quadrupole mass moment 
J2. For a moderate eccentricity its analytical expression is 

where R is the Earth's equatorial mean radius. 

In j26j it is proposed to launch LARES in the so called frozen-perigee 
orbit characterized by the critical value of the inclination, i = 6.34 • 10 1 
deg, for which the ^-precession of (j39j) vanishes. Moreover, it seems that 
Ciufolini suggests to use the LARES data together with those from LAGEOS 
and LAGEOS II in order to measure the Lense-Thirring as well. 

We will now show that such proposals are unfeasible. 

6.2.1 The impact of the even zonal harmonics of the geopotential 

Indeed, apart from the fact that the unavoidable orbital injection errors 
would prevent to exactly insert LARES in orbit with the required incli- 
nation, it turns out that the impact of the other uncancelled precessions 
induced by the even zonal harmonics of higher degree, along with their sec- 
ular variations, would totally swamp the LS effect. In Table El we use the 
calibrated sigmas of the even zonal harmonics of the recently released com- 
bined CHAMP+GRACE+terrestrial gravity EIGEN-CG01C Earth gravity 
model [29 j in order to calculate the mismodelled residual classical precessions 
on the perigee of LARES by assuming i = 6.34 • 10 1 deg. It can be easily 
seen that the mismodelled precessions induced by the first seven even zonal 
harmonics are larger than the LS precession. Due to the extreme smallness 
of such an effect, it is really unlikely that the forthcoming Earth gravity 
models from CHAMP and GRACE will substantially change the situation. 
This rules out the possibility of using only the perigee of LARES. 

6.2.2 The linear combination approach 

In regard to the possibility of suitably combining the Keplerian orbital 
elements of the existing LAGEOS satellites and of the proposed LARES 



~~dl 
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Table 5: Mismodelled classical secular precessions, in mas yr , of the 
perigee of an Earth satellite with semimajor axis a = 1.2270- 10 4 km, inclina- 
tion i = 6.34 • 10 1 deg, eccentricity e = 4.0 • 10 -2 according to the combined 
CHAMP+GRACE+terrestrial gravity EIGEN-CG01C Earth gravity model 
up to degree t = 20 (1-sigma). For the precessions induced by the secular 
variations of the even zonal harmonics, referred to an observational time 
span of one year, the values for Ji,£ = 2,4,6 of j5U] have been used. Such 
precessions grow linearly in time. In the last three rows a - has been in- 
serted because the corresponding mismodelled precessions amount to 10~ 4 
mas yr -1 . The LS secular precession amounts to =f4 • 10 -3 mas yr -1 for an 
object orbiting a central mass along a circular path. 



Even zonal harmonic 


Mismodelled precessions (mas yr 1 ) 


J 2 


1.26 • 10" 1 


■h 


4 • 10~ 3 mas yr -2 


Ja 


4.390 




1.406 mas yr -2 


■h 


1.355 


■k 


6.63 • 10 _1 mas yr -2 


■h 


2.08 • lO" 1 


</io 


2.5 • 10~ 2 


Jl2 


2.5 • 10~ 2 


Jl4 


7- 10~ 3 


Jl6 




Jl8 




J20 
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|31l 1321 I33| in order to reduce the impact of the systematic errors of the 
geopotential on the proposed measurement, it is unfeasible as well. 

Indeed, the perigee is also affected, among other things, by the Einstein 
precession PU, whose nominal value for LARES is 3.280136 • 10 3 mas yr , 
and by the Lense-Thirring effect, which, for i = 6.34 • 10 1 deg, nominally 
amounts to -4.1466 • 10 1 mas yr -1 . This means that if we want to measure 
the LS precession independently of such quite larger Newtonian and post- 
Newtonian effects we would need ten Keplerian orbital elements in order to 
write down a linear system of ten equations in ten unknowns (the first seven 
even zonal harmonics, the LS effect and the two relativistic precessions) 
and solve it for the LS precession. Instead, we would have at our disposal, 
in principle, only four reliable Keplerian orbital elements 5 : the nodes Q 
of LAGEOS, LAGEOS II and LARES-which are affected by the Earth's 
geopotential and by the Lense-Thirring effect but not by the LS force-and 
the perigee of LARES. They would only allow to cancel out the general 
relativistic effects and just one even zonal harmonic. 

The data from an hypothetical satellite with a = 3.6- 10 4 km, as proposed 
in 26 , could not be used for the following reasons. If, on the one hand, the 
classical geopotential precessions would be smaller than the LS effect, apart 
from those induced by J4 and Jq according to EIGEN-CG01C, on the other 
hand, the gravitomagnetic Lense-Thirring and the gravitoelectric Einstein 
post-Newtonian precessions on the perigee would amount to -1.8512 mas 
yr -1 (for i = 6.34 • 10 1 deg) and 2.409956 • 10 2 mas yr -1 , respectively. This 
means that, even for such a higher altitude, one could not analyze only 
the perigee of LARES whose Keplerian orbital elements should be, instead, 
combined with those of the existing LAGEOS and LAGEOS II: with the 
nodes of the LARES, LAGEOS and LAGEOS II and the perigee of LARES 
it would be possible, in principle, to disentangle the LS effect from the post- 
Newtonian precessions and one classical even zonal harmonic. For such a 
high altitude the period of the node of LARES would amount to ~ 10 4 
days, i.e. tens of years. The tesseral K\ tidal perturbation, which is one 
of the most powerful harmonic time-dependent perturbations which are not 
cancelled out by the linear combination approach, has just the period of the 
node. Then, it would act as a superimposed linear bias over an observational 

5 Note that the other routinely and accurately laser-tracked SLR satellites which could, 
in principle, be considered are Ajsai, Starlette and Stella: the useful orbital elements are 
their nodes and the perigee of Starlette. However, since they orbit at much lower altitudes 
than the LAGEOS satellites they would practically be useless. Indeed, they are sensitive 
to much more even zonal harmonics of the geopotential so that they would greatly increase 
the systematic error induced by them. 
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time span of a few years. Moreover, as shown in |35M36| . when high-altitude 
satellites are included in linear combinations involving also lower satellites as 
the existing LAGEOS and LAGEOS II it turns out that the orbital elements 
of the higher SLR targets enter the combination with huge coefficients which 
amplify all the uncancelled orbital perturbations. This would also be the 
case for the K\ tide affecting LARES. 



6.2.3 The impact of the odd zonal harmonics 

The perigee of an Earth artificial satellite is also affected by long-period 
harmonic perturbations induced by the odd {I = 3,5,7...) zonal (m = 0) 
harmonics of the geopotential. The largest perturbation is induced by J3: 
it has a sinusoidal signature with the period of the perigee. Its analytic 
expression is [37| 



dui 
~dt 



-n 



J3 



Jcj sin lo 



esini(l — e 2 ) 3 



+ 



1 



35 



sin 2 i cos 2 i 



• sin 2 i — 1 I sin 2 i+ 



(40) 



For i = 6.34 • 10 1 deg the period of the perigee of LARES, given by (|39|) . 
is of the order of 10 5 days; moreover, the second term of the right-hand- 
side of (|40[) does not vanish. This means that the perigee of LARES in the 
critical inclination would be affected by an additional semisecular bias due 
to J3 which, over an observational time span of some years, would resemble 
a superimposed linear trend. According to EIGEN-CG01C, its mismodelled 
effect would be < 5 • 1CP 1 mas yr _1 . 

This additional bias should be accounted for both in the perigee-only 
scenario and in the linear combinations scenario. Note also that for a = 
3.6 • 10 km and i = 6.34 • 10 1 deg the period of the perigee would amount 
to ~ 10 7 days, i.e. ~ 10 4 years. 



7 Conclusions 

The Dvali-Gabadadze-Porrati braneworld model, in the Lue-Starkman ex- 
tension related to a spherically symmetric central mass, is very interesting 
because it predicts, among other things, small modifications of the Newton- 
Einstein gravity in the weak-field approximation which have testable phe- 
nomenological implications over the Solar-System lengthscale. 
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In this paper we have explicitly worked out its effects on the Keplerian 
orbital elements of the orbit of a test particle without restricting to circular 
orbits. It turns out that the pericentre to and the mean anomaly Ai undergo 
secular precessions which are independent of the size and the shape of the 
orbit. The sizes of these rates are ~ =f4- 1CT 4 " cy _1 and ±1.4- 10~ 3 " cy _1 , 
respectively. The first nonvanishing corrections due to the eccentricity are 
of order 0(e 2 ). They are of the order of 10 -5 " cy _1 . 

The possibility of observing such effects in the orbital motions of the 
inner planets of the Solar-System has been examined, with particular em- 
phasis on Mars. The mean longitude A has been considered along with 
various competing Newtonian and Einsteinian effects. For A the braneworld 
shift amounts for all planets to ~ 1 mas over almost one century while 
the present-day observational accuracy, based on the processing of almost 
300000 data of various kinds for the EPM2004 ephemerides spanning a 92- 
years temporal interval, is 3 • 10~ 3 mas for Mars. A suitable linear combina- 
tion with Venus would allow to reduce the impact of the systematic errors. 
The observational error would be ~ 7%. 

The possibility of measuring the LS pericentre precession in the gravita- 
tional field of the Earth with the perigee of the proposed LARES/WEBER- 
SAT satellite, although appealing, cannot be realized because of the system- 
atic errors due to the mismodelling in the even zonal harmonics of the ter- 
restrial gravitational field. Since the L ARES/WEBER-S AT is conceived to 
be a passive satellite, without any active mechanism of compensation of the 
non-gravitational perturbations, also such kind of systematic errors, which 
the perigee of the geodetic satellites are particularly sensitive to, would be 
fatal. 
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